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Abstract. In the spirit of the generahzed one-particle density matrix for fermions, we 

introduce generahzed one- and two-particle density matrices to state representability condi- 

^ , tions up to second order for boson systems without assuming particle number-conservation. 

lO ' Furthermore, we show for both particle species that, for a semibounded Hamiltonian, the 

^^ \ infimum of the variation of the energy functional w.r.t quasifree states coincides with the 

one of a variation over pure quasifree states. Moreover, it is proven for fermions that only 

pure quasifree states have a generalized 1-pdm that is a projection, and a similar statement 

"^ ' for bosons. 
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1. Introduction 

The Rayleigh-Ritz variational principle for the ground state energy is the starting point of many 
computations and approximations in quantum chemistry. For a many-particle system whose 
dynamics is generated by a Hamiltonian H, it can be written as 



Egs 



inf{tr^(pH) |p>0, tr^(p)=.l}, (1) 



where p varies over the density matrices on the Fock space J-^ = J-^[i)] of the system {Egg 
in ([T]) is actually the total ground state energy in the grand canonical ensemble). A typical 
many-particle Hamiltonian is given as the sum H — Ih + V of the second quantization Ih of a 
one-particle operator h and the second quantization V of a pair potential V. Since Ih is quadratic 
and V is quartic in the field operators, one can rewrite ([T]) in terms of the one-particle density 
matrix jp £ >C+(f)) and the two-particle density matrix Tp S £^(f) (g) f)) of a given density matrix 
p e ^+{J-) as 



Egs = 



mi{£{jp,Tp) \p>0, tr^(p)-l}, 
where the energy functional £ is defined by 

■^(TpiTp) := tit^ih-jp) + ^trtj^t, {hVp). 

The computation of the ground state energy and the corresponding ground state vector of a 
quantum mechanical many-particle system is a complex, if not impossible, task and one resorts 
to approximation methods. The Hartree-Fock approximation is one of the first approximations 



that emerged from ground state computations in quantum chemistry [31 [31 H]. In its original 
formulation, the Rayleigh-Ritz principle for the ground state energy in terms of wave functions, 

^g, = inf{(*,H^)^| *e J-, ||*|b = l} 

of a fermion system with Hamiltonian H is replaced by a variation over Slater determinants. 



^HF = inf|($,H$)^ 



Ne 



$ 



fi A- ■■ AifN, {'fi,'fj)[^^Sij\, 



(2) 



where the Hamiltonian H conserves the particle number, i.e., [H, IN] = 0, with IN being the 
particle number operator. The density matrix p — \^){^\ associated to a Slater determinant is 
a pure, particle number-conserving, quasifree state and ^ can be rewritten as 

_Bfjp = inf < tr jr{pM) p is a pure, particle number-conserving, quasifree density matrix >. 

Since the one-particle density matrix 7 of a fermion Slater determinant $ = (/3i A • • • A (pN is 
the rank-A^ orthogonal projection onto span{(pi, . . . , (Pn} and its two-particle density matrix is 
given as r = (1 — Ex) (7 (g) 7), the Hartree-Fock energy can be written as 

^HF = inf{£(7,(l-Ex)(7®7)) 7 = 7^ = 7*, tr(,(7) < ooj. 

In case of purely repulsive pair potentials V, Lieb's variational principle [TH [31 [S] asserts that 

Ehf = inf{£(7,(l-Ex)(7(g)7)) < 7 < 1, tr„(7) = A^}. 

Going back to a description on the Fock space, Lieb's variational principle reads 

^HF = inf -^ tr_7r(/oII) p is a particle number-conserving, quasifree density matrix > , (3) 

i.e., it asserts that the pureness requirement of the quasifree density matrix can be dropped. 
As was shown in ^ , the property [p, IN] = of particle number conservation is also obsolete 
for repulsive pair potentials V, and the Hartree-Fock energy £'hf agrees with the Bogoliubov- 
Hartree-Fock energy E^bhf defined by 



Ebuf '■= uii <tTjr{pM) p is a quasifree density matrix > . 



(4) 



Our first main result is a generalization of Lieb's variational principle ([3]) in several ways. Namely, 
we show that the infimum in ^ is already obtained from a variation over pure quasifree density 
matrices, 

^BHF = -^bhf'' ^ inf -^ tr_7r(pll) p is a pure, quasifree density matrix >, 

under the mere assumption that H is bounded below. Neither repulsiveness of the pair potential 
V nor the form H = Ih -I- V or even the conservation of the particle number by H is assumed. 
Furthermore, we show that -Ebhf ~ ^bhf'' ^'-"" both fermion and boson systems. The precise 
formulation of this first result and its proof is given in Theorem 14.11 Note that, especially for 
boson systems, it is crucial that our result does not require the Hamiltonian to conserve the 
particle number because for most physically interesting models such an assumption would not 
be fulfilled. 



The above result, i.e., Theorem 14. II brings pure, quasifree density matrices p into focus. These 
are fully characterized by their generalized one-particle density matrix 7p defined in terms of 
their two-point correlation functions as 

(/i®/2, 7p(5i®52)) :- trJp[a*{g,)+a{g^)\ [a{h) + a* Q ^)\) , 

where {a*(/), a(/)| / G f)} are the usual boson or fermion creation operators on F^ fulfilling the 
canonical commutation or anticommutation relations, respectively, on T , with a(/) annihilating 
the vacuum and / i— > J(/) =; / being (a fixed antilinear involution on fj which we refer to as) 
the complex conjugation. Here, we implicitly assume tvjr[pa{f)) = 0, for all / G f), i.e., that p is 
centered. This assumption is irrelevant for fermion systems and made without loss of generality 
for boson systems, as is explained below. The higher correlation of the (centered) quasifree 
density matrix p can be computed from sums over products of the two-point correlation function, 
i.e., in terms of 7p, using Wick's theorem. It is well-known [5l [17] that, as a 2 x 2 matrix with 
operator-valued entries, 7p can be written as 

7p = g li'^J, 7p = 7;, ap = ±<, (5) 

where "-}-" holds for boson and "— " for fermion systems, and A :— JAJ denote the complex 
conjugate and A^ :— A the transpose of a bounded operator A £ B{i)). It is easy to check that 

7p > 0, for boson systems, (6) 

< 7p < 1, for fermion systems. (7) 

We restrict our attention to density matrices with finite particle number expectation, for which 

trf,(7p) = tr^(piN) < oo. 

In this case, it is well-known that the converse of ^ and ([7]) holds true in the sense that, given 
7 as in ([5]), with 7 = 7* e £^(()), 7 > 0, a = ±a^ and obeying 

7 > 0, for boson systems, 

< 7 < 1, for fermion systems, 

there exists a centered quasifree density matrix p G C\_ ( J^) such that 

7 = 7p- 
It is furthermore well-known [SJ[T7] that, if p is a pure, quasifree density matrix, then 

7p = 7 , for fermion systems and (8) 

7p — —^pS^p, for boson systems, (9) 

where 

Our second main result is the converse statement: If the generalized one-particle density matrix 
7p of a density matrix p fulfills (O, in the fermion case, or ©, in the boson case, then the density 



matrix is a pure, quasifree state. Note that the quasifreeness of p is asserted, not assumed. The 
precise formulation of this result is given in Theorem 15.11 

As our third result, we derive representability conditions on the two-particle density matrix 
Tp of a boson density matrix p. Similar to G-, P-, and Q-Conditions for fermion reduced density 
matrices, these conditions follow from the positivity 



tY:r{pP;{a*,a)P2{a*,a)) > 



(10) 



of the density matrix p on positive observables of the form P2 (a* , a) P2 (a* , a) , where P2 {a* , a) 
is a polynomial of degree 2 or smaller in the creation and annihilation operators. A crucial 
difference, however, is that p is not assumed to be particle number-conserving, as this would not 
be a fair assumption for boson systems. Hence, the reduction of the general condition pUj) to 
simpler conditions like G, P, and Q is not as straightforward as in the fermion case and is, in 
fact, not carried out in this paper, but is subject to future work. 



2. Second Quantization and Bogoliubov-Hartree-Fock Theory 

Let (f), (•,•)[,) be a complex separable Hilbert space with the inner product (•, •). : [}x[}— >C. For 
any A^ g IN, the A^-particle Hilbert space representing a physical system of A^ indistinguishable 
particles is given as the A^-fold tensor product of copies of [}, i.e., 



N 



i) 



•S)N 



f). 



fe=l 



N 

The inner product (•,•)„«„ : [)«^ x [)»^ ^ C is given by (/W,g(~))^«„ := Jl (/fc,5fc)f, for 

fc— 1 
any /'^^ = /i (g) • • • ® /jv, g^-^^ = gi (g) • • • (g) gjv e f}**^ and extension by linearity. 
The Fock space J-" is defined as the direct sum of all A^-particle Hilbert spaces. 



^^J-[(,]:=0!,«^. 



Af=0 

Here, by convention 1)®° := C and {f^°\g^°'>)f^m ■= 7^°^5^°^ for /^°',.9^°^ e i)^°'> . Any vector 
\E' e J^ can be written as a sequence of A^-particle wave functions /*^^-' G f)®^: 

V Jn=o 

The vacuum vector fl := (1, 0, 0, ... ) S T, is considered as the basis vector of t)®*^. With the 
inner product (•, ■)jr : J^ x J^ — > C defined by 



(M.,$>^:=5:(/W,5 



JV=0 



JN) 



[,»w 



for any ^ = (/^^^)jv_o ' ^ ~ (d^^^) n-o ^ ■^' ^^"^ Fock space is a Hilbert space. 
The particle number operator is defined as 



JV=0 



For any bounded operator B g B{i}), T{B) :— ®^^qB^^ is an operator on J^. In particular, 
T{B) is trace class, r{B) e C^iT), if B e C^(i}) and, for bosons, additionally ||-B||g([,) < 1. 
A detailed description of the Fock representation can be found in [71 |H1 HI HSl H] • 

2.1. Bosons 

The boson Fock space is the symmetric subspace of the Fock space J^, i.e.. 



J-+^^+[t)]:=5 0[)« 



N=0 

Here, the symmetrization operator S G 13{J-') is defined by 

/ 1 ^ \°° 

with Z*^^) = /i (8) ■ • • ® /tv G f)®^, iV G IN, where ©at denotes the symmetric group with 
permutations tt of iV elements. 

Definition 2.1. For any / G f|, the boson creation and annihilation operators are denoted 
by a*(/) and a{f), respectively. Their domain, which lies dense in J^+, is 'D{¥i^) n J^+ — 

I* = (/(^')^^o G J"+ X;w=o (^ + 1) ll/^^-'ir < ooj. A complete characterization of a* and 
a is given by the properties 

and the canonical commutation relations (CCR) 

[a*{f),a*{g)]=0, [a(/), a(g)] = 0, and 
Hf),a*{g)] = {f,g)t^ 

for any f,g £ f), where Ijr g 'S(J^) is the identity operator on the Fock space and [^,5] := 
AB — BA the commutator. 

The creation operator a*(/) is linear in /, while the annihilation operator a(/) is antilinear. 
Furthermore, the creation and annihilation operators are adjoints of each other, a*{f) = (a(/))*. 
Henceforth, we use the abbreviations a^ = a*{ipk) and a^ = a{ipk) for a fixed, but arbitrary 
orthonormal basis (ONE) {(pk}kLi of ^■ 

Unlike the fermion case, the space generated by all boson creation and annihilation operators 
cannot be used to define a C*-algebra. To this end, we introduce the Weyl operators and 
construct the CCR algebra. For a detailed survey, see, e.g., !9!. 

For every / G (), we define the field operator $(/) : !?($(/)) (- J^+ ^ J^+ by 

<I>(/):=i=(a*(/)+a(/)). 

The field operator is essentially selfadjoint on J^+. Therefore, its closure is selfadjoint and we 
denote it by $(/), as well. 

Definition 2.2. For every / G f), the unitary transformation W(/) : J^+ — ^ J^+, called Weyl 
operator, is defined by 

W(/):=exp(i$(/)). 



The Weyl operators satisfy W(/)* = W(— /) and the Weyl commutation relations 

W(/)W(5) =e-5im(/.9>hW(/ + 5) 

for any f,g G f). The commutator of two Weyl operators is completely determined by the Weyl 
commutation relations. Furthermore, we have W(0) — ljr+ . 

A given field operator $(/) with / S f) is transformed by the Weyl operator W{g), 5 e f), as 

W(5) $(/) Wig)* = $(/) - Im (g, /)„ 1^. 

Hence, W^ = W(i-\/2g) defines a unitary transformation, called Weyl transformation, for any 
g € i)- For any / G f), this transformation yields 

W,a*{f)W;^a*if) + {gJ)^ and Wg a{f)W; ^ a{f) + {f, g)^ . 

Definition 2.3. The C*-algebra W generated by |w(/) / G f]| is called Weyl algebra or CCR 
algebra. 

This algebra is unique up to ^-automorphisms (Cf. Theorem 5.2.8. of [5]). 

Boson Bogoliubov Transformation 

Remark 2.4. The following definition of the Bogoliubov transformation depends on the choice 
of the ONB {ipk}'kLi of ^li since the complex conjugate of a function f £ i), that is given by 

f — Yl i^Wk with some /Xfc G C, fe G IN, is defined by 

k=l 

00 

J ■.= Y^Tlk^k. (12) 

fe=l 

Furthermore, we define for any operator A the complex conjugate operator A by 



{f,Ag):^{f,Ag). (13) 

We emphasize that there is also a formulation to obtain a basis independent definition of the 
Bogoliubov transformation, e.g., in 18, 15_. There, the underlying space is f)® [)* instead of f}0[) 
and an antilinear map J : [)—>■[)*, defined by Jg{f) '■= {g, /)(, for any f,g(z f), and its inverse 
J* : f)* — 7> f) are required. Then, the second component of a vector /®giGf}ffi()is replaced by 
Jg G f)* such that f ® Jg € () ffi t)*. The new vector is antilinear in the second component. This 
supersedes the definition of the complex conjugate of a function. Furthermore, some operators 
map from [}* to f) or vice versa, e.g., v : ()—>■()*. Then, for instance, v and u of the following 
definition are replaced by the maps v : [} — ^ t)* and JuJ*: f}* — )■ f)*, respectively. 

For any linear operator A, the transpose is defined as A^ := A —A*. Now, we are prepared 
to define a boson Bogoliubov transformation. 

Definition 2.5. A linear map U = {^^) : ()©[}— >-i)®i} is called boson Bogoliubov transfor- 
mation if the linear operators u : f) — ?► f) and «:()—?►() fulfill 

uu*-'VV* = lff, u*u—-v ?7=lf|, (14a) 

u*v — ti u = 0, uv — vu =0. (14b) 



Remark 2.6. Any boson Bogoliubov transformation U is invertible. The inverse is given by the 
Bogohubov transformation 

u-^ = su*s 

with 

Remark 2.7. Eqs. (|14b,.b) on u and v are equivalent to stating 

U*SU^S, USU*=S. 

Lemma 2.8. Let U = C^ ^) : f)0f)— >f)ffi() be a boson Bogoliubov transformation. There is a 
unitary transformation Uf/ : T^ — > T^ such that 

^U [a* (/) + a{g)] U^ = a* [uf + vg) + a{vj + ug) 

for all f,g&h'i-f o.'>^d only if v is Hilbert- Schmidt. We call U[/ unitary representation or 
implementation of U on F^ . 

The condition that v is Hilbert-Schmidt is named after Shale and Stinespring |16) . 

States and Density Matrices 

Next, we introduce the notion of states and, afterwards, of density matrices. 

Definition 2.9. A continuous linear functional lu e W* on the CCR algebra W is called a 
state if it is normalized and positive, i.e., w(ljr) = 1 and uj{A) > for all positive semi-definite 
operators A g W. 

Since the boson creation and annihilation operators are not in W, their expectation values 
are not well-defined for all states. In order to find well-defined expressions for these expectation 
values, first, we restrict ourselves to specific states and, then, extend the domain for these states 
appropriately. 

Let W(/) denote a Weyl operator for any / G f} and let oj be a state. We assume that the 
map Tf : R ^- C, t i-^ a;(W(t/)) is four times continuously differentiable for all / S f), shortly 
Tf e C^(]R;C). This assumption provides the definition of the expectation value of a single 
creation or annihilation operator and of the particle number operator. E.g., we have 



^(*(/)):=^^(W(t/)) 



< oo 
t=o 



for any / G f) and, hence, by linearity of uj 



a.(a(/))^-i=[c.(<I>(/))+ic.($(i/))]. 

Analogously, we give a meaning to 

^{a*{f)), a;(e(/)e(g)), and a;(e(/i) e(/2) 6(52) e(.gi)) 

for f,g,fi,f2,9i,92 G f), due to Tf G C"'(R;C). Here, e denotes either the creation operator 
a* or the annihilation operator a. Thus, the expectation value in this state a; can be defined 



not only for elements of the CCR-algebra, but also for polynomials of degree 4 in creation and 
annihilation operators. In order to exemplify such polynomials, we note that, for instance, a 
general polynomial of degree 2 can be written as 

N N 

k,l = l k=l 

with an ONB {ipkj'^^i of f), some iV e IN, and coefficients afe;,^fc;, e^;, Cfci'ffe; M ^ C. 

Definition 2.10. We denote the closure of the indicated extension to the polynomials of degree 
4 in creation and annihilation operators by A'^ . 



For any ONB {y^k}k^i of (), the monotonously increasing sequence of the polynomials 



IN ■ = 



N 



J2 flfcOfc, iV e IN, converges strongly to the particle number operator IN on I?(]N) n T^ 



k=l 



Definition 2.11. Let w be a state. If T/ e ^^(11; C) for any / e f) and a;(]N2) := lim lo{%%) < 
cxD, we write w G Z^ . 

Remark 2.12. For any ui G Z^ , the Cauchy-Schwarz inequality yields 

w(M) < \/uj{%^) < 00. 
Definition 2.13. A state cu E Z+ is called pure if there is a ^ G J^+, such that, for any A G .4+, 

Definition 2.14. A centered state is a state w G Z+ with 

oj{a*{f))=0 (15) 

for any / G t). We denote the set of all centered states by Z^^^. 

As follows from (fT5|) . we also have uj(^a{f)) = for w G 2+j^. 

Definition 2.15. A state ui G Z+ is called quasifree, shortly ui G Zi, if there is a positive 
semi-definite operator /i^ on f) and /cj G f) such that, for every / G f), 

c.(W/)-cxp(-2i(^,/),,-(/,(l^+/i^)/)J. 

The subset of pure quasifree states is denoted by Z^ r. 

Remark 2.16. Zi and Z^ ^ are invariant under Bogoliubov and Weyl transformations, i.e., the 
transform of a (pure) quasifree state is (pure) quasifree, as well. 

Definition 2.17. We say that a state ui G Z+ is coherent, uj G 2^j^, if there is an / G f) such 
that for all A e A'^ 



uj{A) = {n,WfAw}n) 



F' 



Remark 2.18. In particular, wc have 

2+f C Z+ and Z+, C Z^^, C Z+. 

Definition 2.19. We call a positive semi-definite operator p E C^{F^) with trjr+ (/?) = 1 a 
density matrix. 



For a density matrix p £ £^(J^+), the map A^ -^ C,A h^ tr_7r+ (p^Ap^ j defines a state. 

In particular, for every state ut € -E+, there is a density matrix p with trjr+ ip^Ap^ J = i^(^) 

for all A G yl+. Therefore, the notions of pureness, quasifreeness etc can be transferred to the 
corresponding density matrix. 

One- and Two-Particle Density IVIatrices and Representability 

For systems with pair-interactions, the formulation of the variational problem can be reduced by 
the notion of one- and two-particle density matrices. 

Definition 2.20. For any state ui E Z+, the corresponding (boson) one-particle density matrix 
(1-pdm) 7tj : f) — > t} is defined by its matrix elements 

(/,7„g)^:=c.(a*(.g)a(/)) (16) 

for every /, 5 € f) . 

Since any state is positive, we have 

(/, 7c. /)„= c. (a* (.f )«(./)) >0 

for any / € f). Hence, the 1-pdm is a selfadjoint and positive semi-definite operator. Moreover, 
7(^ e C^{i}) due to 

00 00 

tr(,(7„) ^'^{(pk,7c^Vk)f, =cj(^^a^afcj = cj(iN) < 00. 
fe=i fc=i 

Definition 2.21. The (boson) two-particle density matrix (2-pdm) F^^: 1)®!)— >f)(8)f)ofa state 
w G Z+ is defined by 

(/i <» /2, r^ (51 <» 92)) :- ^{a*{g2) a*{gi) a(/i) a(/2)) 

for any fi,f2, 91,92 e I)- 

The 2-pdm is a selfadjoint and positive semi-definite trace class operator since 

tr„®„(F^)=cj(iN2-iN) <^ 
and, for any ip = J2T,i=i l^-kiVk -Sxpi e i) (g) i), pu ^ C, 

00 

(V',r^V')(,,8(, = Yl ^^kl^ij^{ala*laJai) = uj{PP*) > 0, 

i.j,k,l—l 

where P := ^j, ^^-^ pkia*j,a*i . Furthermore, the 2-pdm is symmetric, i.e., Fi^ Ex = ExF^^ = Fj^ for 
all /, g G f) . Here, the exchange operator Ex : ()®[)^l)®f)is the linear map defined by 

Ex(/®5) :=g®/ 

for any /, 17 G (). Summarizing the basic properties of the 1- and 2-pdm, we introduce the notions 
of admissibility and representability. 



10 



Definition 2.22. We call a pair (7,r) of operators on [} x ([) ® f)) admissible if 

(i) r e £^ (f) ® ()) is symmetric, i.e., ExF = FEx ~ F, and selfadjoint, and 

(ii) 7 e £^ ((]) with trf, (7) = a;(lN) is selfadjoint and positive semi-definite. 

Definition 2.23. We say that the pair (7, F) of operators on t) x ([) f)) is representable if there 

is a state a; G Z'^ with 7^^ = 7 and V^ = F. 

Necessary conditions on the pair (7, F) to be representable are called representability conditions. 

In particular, every representable pair (7, F) is admissible. 

Generalized One- and Two-Particle Density Matrices for Bosons 

In [5], a generalized 1-pdm is defined for fermions on the space f) ® [}. We provide a definition 
of the generalized 1-pdm for bosons and, then, further generalize the one- and the two-particle 
density matrices. Here, again, the definitions depend on the choice of the ONB of [), since we 
use complex conjugates of functions, as well as operators as explained in Remark 12.41 We refer 
the reader to |18j for a basis independent formulation. 

Definition 2.24. For any state uj £ Z+, the generalized 1-pdm •^^ is an operator on f) [} 
defined by 

((/i © /2) , 7. {91 © 52)) := ^{ [a*{gi) + a{g,)] [a{h) + a*{7^)] ) (17) 

for /i,/2,gi,ff2 e f). 

Remark 2.25. Defining 

«:: f)-^t), {f,alg):^uj{a*{g)a*{7)), (18) 

we are able to write the generalized 1-pdm as 

7w a^ 

a matrix with operator- valued entries. The 1-pdm 7lj is selfadjoint and a^ symmetric, i.e., 
a^ = a^ for its transpose a^ := a* . 

Lemma 2.26. For any w e Z^ , the generalized 1-pdm 7^^ as defined in (|17p is a positive 
semi-definite operator on f) ® f). In particular, it is selfadjoint. 

Proof. By setting gi — f\ and gi — fi^ the assertion is a consequence of (IT7|) and the positivity 
of the corresponding state. D 

Therefore, the boson 1-pdm 7 is positive semi-definite, too. Unlike the fermion case, the boson 
1-pdm is not bounded above by If,. 

So far, the definitions and statements are well established and can be found for example in 
[15j . in [18] for both particle types and, in a version for fermions, in [5]. 

Lemma 2.27. Let lo e Z^ he a state and 7: f)©[}— >■()©[) iis generalized 1-pdm. For a boson 
Bogoliubov transformation U : [}®f)— ^()®() with unitary representation U;/ : J^"*" — >■ J-""*", define 
utij by UJu{A) := ^(lUc/AlU^) for any A G A'^ . Then, the generalized 1-pdm jjj of the state lou 
is given by 

^u^U*jU. (19) 

Furthermore, 757 = —7 implies "fu Sju = —ju . 



11 



Proof. We consider the matrix elements of ju- For the first assertion, we obtain 

((/i © f2),lu (51 ® 92)) = io {Vu [a*{gi) + a{g,)] U^Uy [a(/i) + a*{7,)] V^) 
= w ([a*(wgi + V52) + a{ug2 + w^i)] 

= (f/(/i®/2),7t^(5i®52)) 

for any /i,/2,5i,ff2 G f). Thus, dH]) holds. 

The second assertion follows from (IT^ and USU* = 5: 




lu S ^u — U* J USU* J U = U* J S J U = —U* J U ^ ~ju, 

which completes the proof. D 

In the following, we give a further generalization of the 1-pdm on the space ^gen := t) ® t) ® C. 

Definition 2.28. For any state uj E Z+, the further generalized 1-pdm j^^ : i^gcn -^ -^gon is 
defined by 

{G,%F) := u{[a*{f,)+a{7,)+^^] [a{g,) + a*(3,) +v]) (20) 

for i^ EE /i ® /2 ® ^ and G = 51 ® 52 ® i^ e i^gcn- 

Remark 2.29. We rewrite the further generalized l-pdm 7;^ as a 3 x 3-matrix: 

(21) 
1 " 

Here, the first moment 6^^ G f) and its dual element 6* G ()* are given by 

{g,b^)^:=u;{a{g)) and b^ ■ g = {b^,g)^ = u{a*{g)) (22) 

for every 5 € f). For the complex conjugate bu of the wave function 6(^ G f), we have {g,buj}, — 

{9,buj)f, ^uj{a{g)) =uj{a*{g)). 

Proposition 2.30. The further generalized 1-pdm is positive sem^i- definite and selfadjoint. 

Proof. The selfadjointness is a direct consequence of ((2T|) . By setting F = G in ((20)) . 7 > 
follows from the positivity of the state to. D 

/ 7 a b \ 

Lemma 2.31. Let j — i a* 1^+7 t \ : [)0()®C— ?►[}©()©€ fee a positive semi-definite trace 

\b' b* 1 / 
class operator with 6 G f), 7 G >C^(fl), and a G £^(f)). Then, there is a unique quasifree state uj 
that has 7 as its further generalized 1-pdm. 

In particular, for any positive semi- definite trace class operator "j — (q* 1, +7) '■ f) ® t) — > () ® t} 
with 7 G C^{\)) and a G £^(()), there is an lj E Z^„f with 7 = 7i^. 

Proof. The second part is a consequence of Theorem 11.4 in [TB] or Theorem 1.6 (i) in [TS]. The 
first part follows from the second part due to the fact that a non-centered state with first moment 
6 G t) is completely characterized by a Weyl operator W^ and the centered state ujq defined by 
uJoiA) := w(W*AWb) for a.ny A £ A^ . D 
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Analogously, we define a generalized 2-particle density matrix F on 
Technically, the generalized 2-pdm should be defined on 

^gen «) -SOgcn =(0f)«)f)]®(0f)]®C. 

It suffices, however, to consider i^sim, since for any polynomial of degree 1 in annihilation and cre- 
ation operators there are an ONB {^k}k=i ^^ t), an iV € IN, and coefficients /i^, i^k, f^kTk, fj-k, I'k G 
C,fc = l,...,iV, such that 

N N 

y^ {^J.kal + Vkttk + akttk + Tfcflfc) = ^ (Afeflfe + t'feOfc) ■ 

fe=l k=l 

For any Af G IN and a given ONB {y^fej^i of f), we set F := (^1,^2,^3, i^4)^ , G := 
(Gi,G2,G3,G4)^, / := {fi,f2f , and 5 := (51,32)^ with F^ := Em=i Mfe/Vfc «"^i, G^ := 
Em=i 4/Vfe «) ^/ € f) «) f), fj ■■= YlkLif^k^'Pk, and ffj := ^^li 4^Vfc e f), where the coeffi- 
cients 4;\ i^^';\ 4-''\ i/^^'^ e C with fc,/ e {!,..., M}, i e {1,2,3,4}, j e {1,2}. Then, we 
define the polynomials Vi and 7^2 by 

M 

k=l 

M 
T> tJ?\ Y^ /^ (1) * * 1 (2) * , (3) * , (4) \ 

V2{F) := 2^ (^/4; Ofeft; + 1^1/ a^^ai + fil/ akUi +fili'akaij . 

k,l = l 

Definition 2.32. The generalized 2-pdni F^ is defined by 

for any F, G G 0^ (f) ® f)) with Y.Z1 f^kk ^J2Zi ^Ifc < °°^ /, g G f) ® f), and a*, ^^ e C as an 
operator on ijsim- The polynomials Vi and 7^2 are of the form specified above. 

As for the generalized 1-pdm, an easy consequence of the definition are the following properties. 

Proposition 2.33. The generalized 2-pdm is selfadjoint and positive semi- definite. 

An explicit form of the generalized 2-pdm as a 7 x 7- matrix is given in Appendix |21 

2.2. Fermions 

The fermion Fock space J-^ = T^ [f)] is defined to be the orthogonal sum 

00 



-^~[fl]:=0[) 



N=Q 
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where, for A^ e IN, 



is the antisymmetric tensor product of A^ copies of () and t) ' :— C Here, the antisymmetrization 



operator A € B{F), A := 0^^^ ^a' with An : f)®^ ^ fl®^, is uniquely defined by 



AnUi®---® In) := 4 E (-1)" h®---®fN=- -4=h A---AfN, 



m ^ ' ' ''^ "'" ^m- 



for /i, . . . , /jv G f), where (—1)'^ denotes the sign of the permutation vr G ©at. 

Definition 2.34. For any / G f), the fermion creation and annihilation operators are denoted 
by c*(/) and c(/), respectively. They are bounded operators on J-~ . Introducing the anti- 
commutator {A, i?} :— AB + BA, they are completely characterized by the properties 

c(/)r!-o, c*(/)r! = /, 

and the canonical anti-commutation relations (CAR) 

{c*(/),c*(g)} = 0, {c(/),c(g)} = 0, and 
M/),c*(g)} = (/,.g)^l^ 

for any f,g e\). 

Definition 2.35. The C*-algebra y4^" generated by {l,c*(/), c(/)| / G f)} is called CAR algebra. 

Let {<Pfc}j.^i be a given ONB of f) and, for this basis, c^. = c*((pfc) and c^ = c{ipk)- For any 
A^ G IN, an ONB of A^-particle Hilbert space f)^^ is given by 



=L^ 



1 < fci < • • • < fc 



AT 



{4. 

Moreover, 

{cl^ • • • cl.„r!| A^ G IN U {0}, 1 < fci < • • • < fcjv} 
is an ONB of the fermion Fock space T~ . 

Fermion Bogoliubov Transformation 

In this section, we fix an (arbitrary) orthonormal basis {"ySfclj-^i of f). As the definitions of 
complex conjugates of both a wave function and an operator depend on the choice of the ONB 
of f), so do the transforms defined in the following. We refer the reader to [18] for a basis 
independent formulation. 

Definition 2.36. A linear map U = {^^) : f)®f)— ?►()©() is called fermion Bogoliubov 
transformation if m : f) — > f) and w : t) — > f) are two linear maps fulfilling 

uu* + vv* = tt,, u*u + v tJ=l(,, (23a) 

u*v + v^u = 0, uv^ + vu^ = 0. (23b) 
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Remark 2.37. Eqs. (j23b .b) on u and v are equivalent to the condition that U is unitary, i.e., 

Therefore, the inverse of a fermion Bogoliubov transformation U exists and is the fermion Bo- 
gohubov transformation U*. 

Lemma 2.38. Let C/=(^S)- f)®f)— >f)®f)&ea fermion Bogoliubov transformation. There is 
a unitary transformation Uf/ : T^ — > J^~ such that 

Vu [c*{f) + c{g)] VI = c*iuf + vg) + civj + ug) 

for all f (B g (z i) ® i) if and only if v is Hilbert-Schmidt. 

This condition on v is called Shale-Stinespring condition [16 . The proof of this lemma can be 
found for instance in jlj- 

States and Density Matrices 

Definition 2.39. For fermions, states are continuous linear functionals w <E {A~)* on the CAR 
algebra which are normalized, uj{\jr) = 1, and positive, lo{A) > for all positive semi-definite 
operators A e A^ . 

Since the fermion systems considered in this work — like atoms and molecules — are particle 
number-conserving, we only deal with even states, i.e., for every odd n e IN, we have 

c^(e(/i)---e(/„)) = 0, 

where e denotes either a creation operator c* or an annihilation operator c. Furthermore, we 
want the particle number expectation value and variance to be finite. Thus, we restrict ourselves 
to the following subset: 

Definition 2.40. We denote the set of all even states with finite particle number variance by 

Z^ = lu e (A^y uj is a state with a;(iN^) < oo and a;(e(/i) • • • e(/„)) = V odd n e f| . 

Again, e denotes either a creation operator c* or an annihilation operator c. 
Remark 2.41. By the Cauchy-Schwarz inequality, we have for any state lo G Z^ 

a;(]N) < yC(i^< oo. 
Definition 2.42. A state oj G Z^ is called pure if there is a ^ G J^~ , such that 

uj{A) = (*,A*)^ 
for any A E A^ . 

Definition 2.43. A state lj E Z^ is called quasifree, shortly ui E ZZ, if it fulfills Wick's 
Theorem, i.e., 

w(ei 62 • • • e2Ar-i) = and 

w(ei 62 • • • 62Ar) = ^(- l)''a;(6^(i) e^(2)) • • • uj{e^^2N-i) e^(2N)) (24) 
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for every N £ M, where e^ denotes either a creation or an annihilation operator for every 
i € {1,2,..., 2iV}. The sum is taken over ah permutations tt £ &2N satisfying 

7r(l) < 7r(3) < • • • < 7r(27V - 1) and 7r(2fc - 1) < 7r(2A:) 

for every k £ {1,2, ... , N}. The right hand side of (P^ is cahed Pfaffian. 

The subset of the pure quasifree states is denoted by Z~ r. For any iV G M and any orthonormal 
vectors ipi, . . . , ippf £ f), the vector ipi A ■ ■ ■ A (pN £ J^^ is cahed Slater determinant and defines 
a pure quasifree state. 

Remark 2.44. The sets Z^^ and -Z ~ f are invariant under Bogoliubov transformations. 

There is a characterization of pure quasifree states using the Bogoliubov transformation. 

Remark 2.45. A state uj £ Z^ is pure quasifree if and only if there is a fermion Bogoliubov 
transformation C/:f)®f)— >f)©f) with unitary representation TU : T^ -^ T^ , such that for any 

A£A' 

w(A) = (Ufl,AT[Jfl)_^. 

Remark 2.46. Since we assume that the fermion states are even, they are, in particular, centered 
(see Definition 12.141 for bosons). Moreover, for bosons, the set of centered quasifree states is a 
proper subset of the set of quasifree states (Definition 12. 15L Z^^ C Zi, while for fermions all 
quasifree states are centered. 

Definition 2.47. A selfadjoint, positive semi-dcfinitc trace class operator p £ C^{F^) of unit 
trace, tr^r- (p) = 1, is called density matrix. 



The map A — > C, A !—> trjr- ip'^Ap^ j defines a state. Since we only study fermion systems 

that preserve the particle number, we restrict our attention to density matrices which commute 
with the particle number operator and have a finite squared particle number expectation value, 

oo 

/» = P^'^^ and tr^- (pi %^pi) <oo. (25) 

N=0 

Note that, ii m,n > 0,m ^ n, then 

tr^- {p^c*(/i)---c*(/™)c(5i)---c(.g„)p^| =0 



for any choice of /i, . . . , /„, 51, . . . , 5„ G (), due to (I25j). 

Remark 2.48. In particular, for every state uj £ Z^ , there is a density matrix p fulfilling (j25p 

and ti'jr- (pi Apij ^ uj{A) for aU A £ A^ . 

One- and Two-Particle Density IVIatrices 

We now introduce the notion of fermion one- and two-particle density matrices. 

Definition 2.49. For any ui £ Z^ , the one-particle density matrix (1-pdni) j^ £ /S(f)) of oj is 
defined by 

(/,7-ff)„:=^(c*(5)c(/)) 
for f,g £i}. 
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Definition 2.50. The two-particle density matrix (2-pdni) F^: i) (E) i) —?' i) <E) i) oi a, state a; e Z 
is the bounded operator given by 

(/i ® /2,r^ (51 ® 52))[,55„ := ^(c*(52) c*(5i) c(/i) c(/2)) 

for fi,f2,9i,92 e fl- 

An outhne of basic properties of the fermion 1- and 2-pdni can be found in Lemma 2.1 of [4]. 

Generalized One-Particle Density Matrix for Fermions 

Analogously to the boson case, we define a generalization of the 1-pdm for fermions as in [5]. 
The complex conjugates of a function or of an operator are defined as for bosons in (|12l) and 
(|T^ . respectively. 

Definition 2.51. Let to e Z^ and fix an ONB {^k}^=i of f). Then, the generalized 1-pdm 7^^ 
of oj is an operator on f) ® [) defined by 

((/i © /2) ,7- (gi ® ff2)) := u;{[c*igi) + cig^)] [c(/i) + c*{7^)]) 
for fi,f2,9i,92 e [}. 
Remark 2.52. Again, we define the operator a* : [) — > f) for every /, 5 S f) by 

(/,a:.g):=c.(c*(g)c*(7)). 
Then, the generahzed 1-pdm is expressed as the matrix 



As for bosons, the fermion 1-pdm 7 is selfadjoint, but a is anti-symmetric, i.e., a^ — —a, as 
follows from CAR. 

Lemma 2.53. For any uj G Z+, the generalized 1-pdm 7^^ is a positive semi-definite operator 
on f) © f). In particular, it is selfadjoint. Furthermore, it is bounded above by If, If,. 

We refer the reader to [5] for a proof. From Lemma [2?53] we deduce < 7 < Ip, for the 1-pdm. 
A consequence of Wick's theorem is the following lemma. 

Lemma 2.54. A quasifree state uj G Z^ is uniquely determined by its generalized 1-pdm ^^. 

Furthermore, the generalized 1-pdm transforms in a specific manner under the Bogoliubov 
transformation. 

Lemma 2.55. Let uj G Z^ be a state with generalized 1-pdm 7: f)®f)— 7>f)®f). For a fermion 
Bogoliubov transformation U : [}®f)— >'f)®() with unitary representation Uf/ : T^ — > J-^ define 
utu by UJu{A) :— uj(Vu AJ]^^ for any A G A^ . The generalized 1-pdm ^jj corresponding to the 
state ujjj is given by 

ju = U*jU. (26) 

In particular, 7^=7 implies 7^ —Ju- 
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Proof. For any /i, /2,5i; 32 e (), we have 

((/i ® f2),lu (ffi © 52)) - u;{Vu [c*{gi) + c{M VhVu [c{h) + c*(f,)] U^) 

= ^{[c*{ugi + vg2) + c{ug2 + vg^)] [c{ufi + u/2) + c*{uf2 + vj^)]) 
^{U{gi®g2).lU{h®h)) 

for the matrix elements of 7^/. Thus, (051) holds. Furthermore, by the unitarity of U and (051) . 
we obtain ^Ij = U* j U U* j U ^ U* j'^ U and 7^ = 7 yields jfj = U*^U ^ ju- □ 

2.3. Bogoliubov-Hartree-Fock Theory 

Boson Bogoliubov-Hartree-Fock Theory 

For bosons, the number of particles in most physically relevant models is not fixed. As, for 
instance, in a system of photons interacting with an electron, photons can appear or disappear, 
depending on what is energetically favorable. Thus, the particle number should not be fixed in 
the variational process yielding the ground state energy Egs := inf {cr(H)}. By the Rayleigh-Ritz 
principle, the ground state energy (as well as the ground state) is determined by 

Egs = miiuj{M) uj eZ+\ . 

In the Bogoliubov-Hartree-Fock (BHF) theory, the variation is restricted to quasifree states: 

Ebhf -.^ mi \^Lj{M)\uj e Z+Y 

The BHF energy E'bhf is an upper bound to the ground state energy Egg- Note that, unlike 
the common definitions of quasifreeness, our quasifree states are not necessarily centered. Since 
a quasifree state is uniquely determined by its further generalized 1-pdm 7^^ , there is an energy 
functional Sbhf ■ 2?('?bhf) ^ <C, I?(£bhf) C S([) ® () © C), such that w(H) = £bhf(7c^)- Thus, 
the BHF energy is rewritten as 

^BHF = inf |£'bhf(7cu) w e Z^A = inf |£'bhf(7) 7 > 0, tr(7) < ooj . 

The second equality is a consequence of two facts: On the one hand, any quasifree state uj with 
first moment b is linked to a unique centered quasifree state via the Weyl transformation W;,. On 
the other hand, any positive semi-definite operator j — (a* i(,+7) on f) () fulfilling tr(7) < cx) 
is the generalized 1-pdm of a centered quasifree state, cf. [15) . 

Fermion Bogoliubov-Hartree-Fock Theory 

For fermions, assume U : R'^ — 7> R to be an external potential and V : R'^ x E,'^ — > ]R,q a 
repulsive interaction between two particles. There are multiplication operators associated to 
these potentials which we also denote by U and V, respectively. With the Laplace operator A, 
the Hamiltonian of the system is given by 

^ 1 

i=\ l<i,j<N 
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where a;, G R"^, 1 < i < A^. We only allow for potentials for which Hm is defined as a selfad- 
joint operator on a dense domain Djq and is bounded below. The second quantization of this 
Haniiltonian is 



H 



7 , ^ij "'i '-J 



jCa C-i 



_j oo 



2 zl^ /jj,fc;CjqcfcQ, 



where the one-particle operator h and the interaction operator are given by 

respectively, for any elements of a given ONB {</3i}i=i of f) with || V((9i|L < oo. The Hamiltonian 

Hn is the restriction of H to the A^-particle Fock space /\ f). If we do not assume the dynamics 
to conserve the particle number, the ground state energy of the TV-particle system is determined 
by the Rayleigh-Ritz principle: 

we z-\ . 



Egs = inf |cj(H) 



Using the energy functional 



£ij,T):^tT{hj) + ltT{VT) 



this can be re-expressed as 



Egs = inf < £(7, r) (7, r) is representable > 



Here, the problem of representability arises, i.e., a classification of all representable pairs in 
() X () ® (). In order to obtain an upper bound to Egs, the variation is restricted to quasifree states 
which yields the Bogoliubov-Hartree-Fock energy 



E 



BHF • — 



inf |cj(H) 



uj e Zq(.| = inf |£bhf(7) 7 > 0, tr,,(7) < ooj 



For any quasifree state a;, the Bogoliubov-Hartree-Fock functional £bhf is given by fBHF(7Lj) '■— 
w(H), where j^j is the generalized 1-pdm of uj. 

3. Bosonic Representability Conditions and the Generalized 
Two-Particle Density Matrix 

3.1. Particle Number-Conserving Systems 

To our knowledge, sets of representability conditions given in the literature are for particle 
number-conserving systems for fermions, as well as for bosons. I.e. only states, that fulfill 



CO 



l[a*{h 



.fe=i 



H'^^-^'' 



.1=1 







for any two sets {fk}k=i 1 {ff/}/=i ^ h with tti, n e W U {0} and m ^ n, are considered. 

Since the dynamics of many realistic physical boson systems do not conserve the particle 
number, an alternative should be found. First, we restate some representability conditions for 
bosons. 
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Definition 3.1. Let (7, F) be a pair of operators on [} x (t) (8) ())■ We say that (7, F) satisfies the 
representabihty conditions up to second order with particle number-conservation if 

1. (7,F) is admissible, 

2. F satisfies the P-condition, i.e., 

F>0, 
and 

3. the G-condition, i.e., for any A E S([}) we have 

tr {{A* (g) A) [F + Ex (7 (g) I,,)]) > |tr (A 7)]^ . 



These conditions can be found, e.g., in [T2J1I3]. Note that these conditions are only necessary 
conditions, but do not ensure that the considered operators are one- and two-particle density 
matrices. Furthermore, we omit here other known conditions like the Ti- and r2-condition, cf. 

Remark 3.2. The Q-condition is omitted since it follows from the P-condition and the positivity 
of 7, see [13]. Nevertheless, in the same manner, we can rephrase the Q-condition from [T3] as 

r > - (li,®r, + Ex) (7 ® If, -f If, «) 7 -t- If, (g) If,) . 

The representabihty conditions for bosons up to second order are derived in the same spirit as 
it is done for fermions in |4j. 

Theorem 3.3. Let uj he a linear continuous functional on A'^ such that w(l) = 1, wfW^) < 00, 
and a;(ei . . . e2Ar_i) — for all A^ e IN, where Ck denotes either a creation or annihilation 
operator. Furthermore, let F^j and ^^ he the corresponding 1- and 2-pdm and {fk}'kLi o,''t- ONB 
of f) . Then the following statements are equivalent: 

(i) For any polynomial Vr G A'^ in creation and annihilation operators of degree r < 2, we 
have 

(i'i) 7tj > and F^j fulfills the G- and P-condition. 

Since the proof is analogous to the fermion case considered in [i], we omit the details here. 
Note that, unlike the fermion case, the trace class conditions on the 1- and 2-pdm cannot be 
derived from the polynomials since the boson creation and annihilation operators are unbounded. 

3.2. Systems without Particle Number-Conservation and the Generalized 
Two- Particle Density Matrix 

We generalize the definition of the representabihty conditions up to second order to systems — 
and, thus, states — which do not conserve the particle number. These representabihty conditions 
arise in the same manner as those for particle conserving states by considering expectation values 
of polynomials up to second order in the creation and annihilation operators. Due to the absence 
of particle number-conservation, the expectation values of terms, in which the number of creation 
operators is not equal to the number of annihilation operators, do, in general, not vanish. A 
simple consequence of Definition 12.321 is the following proposition. 
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Proposition 3.4. The representability conditions up to second order are satisfied if the pair 
(7, r) of operators on f) x (f) (g) f)) is admissible and T is positive semi-definite as an operator on 

Let w be a linear functional on the operators on J^"*". Since, on the one hand, any polynomial 
up to second order in creation and annihilation operators can be written asV = 'P2{F) +'Pi{f) + 1^ 
with F e 0"* ()®^ / e 0^ f), i^ e C, Definition 1^3^ yields 




On the other hand, every element of i^sim can be written as a vector with i^ € f)®^, f £ 
[), i^ G C Thus, the representability conditions up to second order are exactly those arising 
from 

uj{VV*) > 

for any polynomial V in creation and annihilation operators of degree r < 2. 

Remark 3.5. Since the generalized 1-pdm appears as a block in the generalized 2-pdm, it inherits 
the definiteness property from the generalized 2-pdm. 

If one varies only over particle number-conserving states, then F assumes a block-diagonal form, 
and the complexity of the representability reduces considerably. In fact, only three independent 
conditions remain which are reminiscent of the G- and P-condition in quantum chemistry (see 
Theorem [ 



4. Variation over Pure Quasifree States and 
Bogoliubov-Hartree-Fock Energy 

For bosons. Theorem 1.2 of [3] states for the Pauli-Fierz model that the Bogoliubov-Hartree- 
Fock energy coincides with the infimimum of the energy functional for a variation over pure 
quasifree states. We prove a more general statement, which holds for bosons, as well as for 
fermions. The main result of this section is the following 



Theorem 4.1. Assume the Hamiltonian H to be bounded below. Then, 
Ebhf = inf <, Lu (H) oj is pure and quasifree > 



ppuro 
-^BHF- 



We show the statement in the following two subsections for bosons and fermions separately. 

4.1. Bosons 

For bosons, a more precise statement of Theorem 14. II is: 

Theorem 4.2. Let H &e a Hamiltonian on T^ that is bounded below. Then, 

FBHF=inf{w(H) 



-eZ+,}^:i?rHF- 



In order to prove the theorem, we need some properties of quasifree and pure quasifree states. 
To this end, we give a characterization of quasifree and pure quasifree states using the Bogoliubov 
transformation. 
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Lemma 4.3. For any quasifree density matrix, there are a positive semi-definite operator C € 
£^([)) with ||C||g/fjN < 1 and second quantization T{C) := ^'^^qC^^ , a boson Bogoliubov 
transformation with unitary implementation HJ, and / G f) such that 

//, furthermore, the density matrix is pure, it is of the form W/ U |rj) (J7| TU* W^, where we used 
the Dirac bra-ket notation. 

This Lemma is a consequence of Lemma IILl in [3]. 



Proof of Theorem \4-S\ Without loss of generality, we assume the Hamiltonian to be positive 
semi-definite. If H is bounded below, there is a constant /i > such that Hq := H + /J,ljr > 0. 
Considering Hq instead of H just adds the constant /i to both E'bhf and -EgHF- 
The inequality 

^BHF = inf {c^(H) I u e Z+} < inf {^(H) | u G Z+ j } = E^"^^ (27) 

follows from the definition of the BHF energy, since the variation is restricted to the proper 
subset Z+qf C Z+. 

It remains to prove that w(H) > -EgHF ^'^^ ^^J quasifree state w G Zi. Let uj £ Zi with 
a;(H) < oo and denote the corresponding density matrix by p. Then, 

w(H) = tr^+ (p^Mp^) = trjr+ ((pi ui) (ui pi\\ 

since H > 0. Therefore, p^ Ha is Hilbert-Schmidt and we obtain by the cyclicity of the trace 
tr^+ ((pi Ui) (ui pi)) = tr^+ (h^ pUi^ 



Since H is selfadjoint, there is an ONB {*fc}^i of J^+, such that *fc e X'(H) for any k eK. 
Then, 



tr^+ (h^pH^) =^(^H5*fc,pH^*fc) . 



fe=i 



By Lemma 14.31 the positive semi-definite operator p can be written as p — kk* , where 



K:=WfV 



{trMncw 



with some / G f), a Bogoliubov transformation with unitary implementation UJ, and some C G 
CHi)),C>0,\\C\\g^^^<l. Hence, 

CXD „ 

w(H) =^|U*H^*fc . (28) 

fe=i ^ 

We continue by introducing a resolution of the identity with coherent states. To this end, we 
consider an increasing sequence of n-dimensional Hilbert spaces f)„ C l)n+i ^ f), n G IN with 
Umpin''" = f) and Cl)n ^ fin- For any n-dimensional Hilbert space [}„, there is an isometric 
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isomorphism /:()„—)■ C". We define the measure d/i„(z(")) on f)„ by J, d/i„(z*^"))/(z'-"') :— 

/c- '^"^T^ f{lz'^''^), where a; := Re(z), y := Im(z). For any n e M, we have () = f)„ e fj;^, 
where ^^ denotes the orthogonal complement of f)„ in [). Moreover, J^+ = J'+[f)„] (g) J"+[();^] 
and fi = ri„ (K) Sl;|; with f2„ S J^+[t)„] and fJ^ e ^"+[1)^]. For every n e M, the projections 
|W(z("))r2) (W(z("))f7|, z(") e ()„, satisfy 






(")^ 



W(z(")) fi) (W(z(")) 17 



see, e.g., [SlII^. Consequently, 

(*,*)^= lim / d^„(z(")) /*,W(z("))r2\ 



ll„ 



J^ 



for any ^ e J^+. Thus, each summand of the right hand side of (j28l) is rewritten as 



K*H5^S 



= lim / d^„(z(")) /H^*fc,«;W(z("))f7\ 



The sequence (fci-^/. d/i„(z(")) (n^^fe, k W(2(")) 17 



,, ^p-n y'' y , jii- I K, •" ■■ V- ; ""/ ; i^ monotonously increasing. 

" ... /n=l 

Therefore, the summation and the limit is exchanged by the monotone convergence theorem, 
where the summation is considered as an integral with the counting measure. Thus, we get 



bj\ 



fe=i"''' 
Afterwards, Fubini's Theorem yields 



OO r, 

H) = lim ^ / dAi„(z("^) /H5*fc,KW(z("))fi\ 



fe=i "^ 

= lim / d^„(z("))(KW(z("))f7,HKW(z("))r2 



J^ 



and we conclude from the proof of Lemma III. 7 in 3 that 



KW(z("))fi = W/U 



(tr^+(r(c))) 



Y W(z(")) f} = ;/c(z*"0 Wg U 17 



for some 5 G t) and vc{z<^^) £ C with lim„^oo J,,^ d^„(z(")) |i/c(2:'"^)| = 1. By Lemma l4?3l 
this vector defines a pure quasifree state and, consequently. 



w(H) = lim / d//„(2(")) (KW(z("))f7,HKW(z("))l7 



T 



> -E^BHF 1™ / dflniz 



,(")^ 



tin. 



vc \z 



M) 



ppurc 
-^BHFi 



which completes the proof. 



D 
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4.2. Fermions 

For fermions, a similar result to Theorem 14.21 holds : 

Theorem 4.4. Let H &e a Hamiltonian on J-^ that is bounded below. Then, 

Ebhf = inf |w(H) lo e Z^qA =: -BgnF- 
Before we prove this theorem, we need two preparatory lemmas. 

Lemma 4.5. Let uj G Z~f with density matrix p. Then, there are a decomposition f) = f)s ® [)r 
with n := dim(f)s) < oo, a positive semi-definite trace class operator B e B{\)t), and a fermion 
Bogoliubov transformation U with unitary implementation TU, such that 

P^^u Ia\vi a • • • a (^„) (^1 A • • • A (p„|) ® ^^^^(pL)J ^*u (29) 

for any ONE {^k}1^i ofi)s. 

Proof. It is known that there are fermion Bogoliubov transformations U such that the generalized 
1-pdm oi pu := V^ pVjj is of the form 

for some < 7c/ < 1(, with tr^^{'-fu) < oo, see, e.g., [SI Theorem 2.3]. Let ()s be the eigenspace 
of 7(7 associated to the eigenvalue 1 with dimension n < oo and f)r its orthogonal complement. 
Then, jjj = Ps + 7r, where Ps is the orthogonal projection on f)5 and 7r the restriction of ju 
to f)r- Note that 7r satisfies f)5 C ker(7r), 7r [)r C ()r, and < 7r < /il() for some < /i < 1. 
Let ipi, . . . ,ipn he an ONB of ()5. Moreover, let 

T{B) 



p' := \ipi A ■ ■ ■ A ipn) {ipiA---Aipn\' 



tr^- (r(i?)) 



with _B := (7r) (Ifjr — 7r) • In order to show that p' = pu, it is sufficient to observe that 
p' defines a quasifree state lo' and j^> — ju from pop . since quasifree states are characterized 
by their generalized 1-pdm, see [5]- Note that we implicitly used the decomposition J-^ = 

Remark 4.6. For any positive semi-definite trace class operator B e 'S(f)r), there are an ONB 
{4>k}kLi of[)r and coefficients 6fc > 0, fc e IN, such that i? = X^fcli ^fe I'f'k) {(pkl and ^^-^ 5^ < oo. 
Thus, 

COO \ OO oo 

(8) nh) = n t"-^- (r(6fc)) = n (1 + ^fc) 
fe=l / fc=l fc=l 

which converges due to X)fc=i ^k < oo- Here, T{bk) should be understood as the second quantized 
operator r(6fe \(j)k} {(f)k\) on J""[C(/)/c]. 

Lemma 4.7. Lei w e Z^ wit/i density matrix p. Then, there is a sequence {pk)k=i '^f pure 



K 



quasifree density matrices and {^k)k=i ^ [0, oo) with X]fc=i ^k < oo such that 

(*i,p*2)^= lim (*i,VAfcPfe*2) 

\ fc=l / ;c- 

for any ^I'l, 4*2 G J^^ . I.e., every quasifree state is a convex combination of pure quasifree states. 
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Proof. From Leninia l4.5[ we know that every quasifree density matrix is of the form (|29p and we 
use the notation specified there in the following. We complete {<fk}^=i to an ONB {<Pfe}feLi of 
f), where {fkj^^n+i i^ ^^ ONB of f)r- Then, 



(*'*)^=1™ 1™^ J2 (*,<^»i A---A¥',:,)^(¥'ii A---A^,;,,$)^ (31) 

fe=0 l<ii<---<ifc<M 

for any ^,$ g J^^ . Choosing the ONB {fk}'kLn+i ^^ 'ir such that B is diagonalized and using 
ni . we obtain in the weak sense 

2 A ^, , r(B^) V 



[ti>- inBw 



N 
M,N^oo 



= .MV" Yl Yl i\ipi A ■ ■ ■ A ifin) {'Pl A ■ ■ ■ A ipn\) 

k=n+l 'n+l<ii<---<ik<M 

T(Bi) T(Bi) 

- \ip^, A ■ ■ ■ A if,,) {if,, A ■ ■ ■ A ip,,\ 



V[tr^-(r(B))]^ [tr^- (^(B))]^ 

This can be written as 

AT 



K' 



" tr^- ivm M^oo Y E (|^iA...A^„)(^iA---A^„|) 



iB^^fi^ A ■■■ AB^^ipiA (B^^ifi^ A ■■■ AB^if,^ 



N 



lim > > 



\.Tt- {V(B)) M,N^oc ^-^ 

■r \ \ 1 1 A;=n+1 n+l<ii<---<ifc<M 

<y5i A • • • A ¥?„ A 55(^,1 A • • • A B^ip.A (ipi A ■ ■ ■ A ipn A Biip,^ A • • • A B^ip 



Each operator 



(pi A ■■■ A(pnAB^ip>i^ A ■■■ AB^ip^A (ipi A ■ ■ ■ A ip>n A B^ip.^ A ■■■ AB^ip,^ is 

either equal to zero or a pure quasifree density matrix (up to a normalization constant). Finally, 
a pure quasifree density matrix conjugated by a Bogoliubov transformation is a pure quasifree 
state, too, which completes the proof. D 

Now, we are prepared to prove Theorem 14.41 Since the proof is, to a large extend, similar to 
the proof of Theorem l4.21 we only give details where there are differences. 

Proof of Theorem \4-4\ Again, without loss of generality, we assume that the Hamiltonian is 
positive semi-definite. 

As for bosons, the inequality 



Ebhf = inf {a;(H) lo e Z^j } < inf {cj(H) 



^e^pqt} =£^BHF 



is immediate. 

Thus, we show w(lE[) > E^p for any to € Z~f. Let lo G Z^ with a;(H) < oo and denote 
the corresponding density matrix by p. Furthermore, let {^fc}^i be an ONB of J^^ , such that 
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^k ^ X>(H) for any A: G IN. Analogously to the boson case, we obtain 

CXJ 



fe=i 



T 



By Lemma 14.51 the positive semi-definite operator p can be written as p = k k* , where 



k:=V 



u 



\ipi A ■ ■ ■ A ipn) {(pi A ■ ■ ■ Aipn\ ' 



[tr^-(r(i?))]^ 



with a decomposition f) = fig ©-'- [)r, n :— dim{l}s) < oo, an ONB {(pk}^=i of t]s, a unitarily 
implementable Bogohubov transformation U, and a positive semi-definite trace class operator 
B e 6(t)r). Hence, 



(H) =^L*H^*j 



fc=i 



Instead of a resolution of the identity by coherent states for bosons, we use the resolution of 
the identity by Slater determinants. 



N 

1^- = lim lim V V \ip,^ A ■ ■ ■ A ipt,) {ipi, A ■ ■ ■ A (p^,\ 

JV— s-oo M-i-oo ^ — ' ^ — ' 

;=0 l<ii<---<ii<M 

as in the proof of Lemma |4. 71 in particular, Eq. pip . Then, we obtain 

N 



fe=l l=Q l<ii<-<ii<M 

(\ oo 
fc i-> limA/_5.oo X) X) (E2 vj/j^^, K((pi^ A • • • A (/3i, )\ is 

/=0 1<ii<---<i,<Af \ '^ J N=l 

monotonously increasing, the monotone convergence theorem allows for a exchange of the k- 
summation and the first limit. Using the monotone convergence theorem a second time to 
exchange the second limit and the fc-summation, we obtain 



oo N 



^W=^1™^J™^EE E |(H^*fc,'*(^nA---A(^j)^ 

fc=i /=o i<n<---<4i<J\/ 

Furthermore, we change the order of the summations, since the sum is absolutely convergent, 
and get 

N 

uj(ll) = lim lim S^ S^ {k (ip^^ A ■ ■ ■ A ipn) ,11 k (ipi^ A ■ ■ ■ A (pn)) ^r ■ 

1=0 l<ii<---<i,<AI 

Every vector k (ipi-^ A • • • A (/Si, ) defines a pure quasifree state, cf. the proof of Lemma 14.71 Since 

(vE.,Hv|/)^><- 
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for any pure quasifree state ^ e J'-, we finally have 



w(H) > ^Pyp lim lim V V (.^,, A • • • A .^,, , p (.^,, A • • • A .^,, ))^ 

i=0 l<ii<---<i!<M 



_ ppurc 
— ^BHF- 

This proves the assertion. D 

Theorem 14.11 now follows from the Theorems 14.21 and 14.41 

5. Pure Quasifree States and their Generalized One-Particle 
Density Matrix 

For a given generalized fermion 1-pdm 7, it is known that there is a pure quasifree state a; which 
has 7 as its generalized 1-pdm, if and only if the generalized 1-pdm is a projection, i.e., 7^ = 7 
(see Sect. 15. 2| ). For bosons, a similar statement is also known. In this section, we show that an 
even stronger relation holds: 

Theorem 5.1. The following statements are equivalent: 

(i) Lo is a centered pure quasifree state. 

(ii) The corresponding generalized 1-pdm 7 satisfies tr (7) < 00 and 

757 = —7 for bosons, 
1 — 7 for fermions. 

Recall S — Ifj Q) (— IrO G 'B(f) ffl ())■ A proof of Theorem 15. II in the boson case is given in the 
following subsection. Two consequences of this theorem are discussed afterwards. In the second 
subsection, we prove the statement for fermions. 

5.1. Bosons 

Before we show Theorem 15. II for bosons, we give some preparatory lemmas. 

Lemma 5.2. // an operator 7 = ( J* 1+^) '■ f) © f) — > f) © f) satisfies 7 > 0, tr(7) < 00, and 

757 = -7, (32) 

there is a centered pure quasifree state uj G -ZJ" £ n 2+„ that has 7 as its generalized one-particle 
density matrix. Furthermore, let lu G ^^qf ^ ^ten ^^ ^ centered pure quasifree state. Then, the 
corresponding generalized 1-pdm 7 fulfills p2p . 

For a proof see e.g. [13 [H]- 

Eq. (|5^ is rewritten in a single equation for operators on [), i.e., we do not need the matrices 
7 and S. 

Proposition 5.3. Let ui G Z^ he a state with the generalized 1-pdm J — {a* 1, +~)- Then the 
following statements are equivalent: 
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(i) ryS-f = --f. 

(ii) 7^ + 7 = a a* . 
Proof. Computing and simplfying the matrix products of (i), we obtain the four equations 

7^ +7 = aa*, (33) 

7^ + 7 = a* a, (34) 

7a = a7, (35) 

a* 7 = 7a*. (36) 

Thus, the imphcation (i) => (ii) is immediate. It remains to prove (ii) => (i). Eq. (j36|) is 
the adjoint of ([55)) . and (IMl) is equivalent to ([55)) . The system of equations reduces to ([55)1 
and (j35p . Furthermore, we show that (|35|) follows from (j33l) . We define / : R+ ^ ]R+ by 
f{y) := -^/y + 1/4 — 1/2 and observe that / is the inverse map of x 1— >■ x + a;^ , R^ ^- R^ . Then 



7 = /(a a*) and 7 = /(a a*) = /(Q!*a). 

Since a a* is bounded, we approximate the function / by a sequence {pn)n=i ^^ polynomials, 

C R. So, 



I"" Hop 'II"" Hop 



i.e., liTan^ooPn{x) — f{x) uniformly on the compact interval 

p„(aa*) and p„(a* a) are well-defined. Using (a a*) a — a {a* a) for all tti € IN and limits 
in operator norm, we obtain 

^a = f{aa*)a= limp„(aa*)a= Ivtn apn{a*a) — a f {a*a) = a"^, 

n— >oo Ti^oo 

which proves the assertion. D 

Remark 5.4. In [3] 

7 = - (cosh(2r) — 1) , a = - sinh(2r) 

are used, where {f,ag)f^ = {f,ag)f^ and r : f) ^ [) is an antilinear operator, r obeys {f,rg) — 
{g,r f) for any f,g & i) and r^ is trace class. These two equations are, however, implied by (j33p 
and, in turn, yield (l35t . 

Centered pure quasifree states can be characterized by a Bogoliubov transformation (see |15] 
for the proof): 

Lemma 5.5. A centered boson state lo € Z^^^ is pure quasifree if and only if there is a boson 
Bogoliubov transformation U : f)®f)— ^-flCDf) with unitary representation Vu : T^ — > T^ , such 
that for any A e A^ 

uj{A) = {Vu^,AVu^). 

The relation between a generalized 1-pdm fulfilling (1321) and the corresponding centered pure 
quasifree state is even closer. 

Lemma 5.6. Let tu G -Z+„ and assume that the corresponding generalized 1-pdm 7 satisfies 

7^7 = -7. (37) 

Then, uj is a centered pure quasifree state. 
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Proof. As stated in Remark 12.251 the generalized l-pdm is of the form 7 = („* 1+7)1 where 
7 : [}->() is the l-pdm and a* : f) — )■ () is defined in (IT51) . For any fc € IN, let 0fc denote an 
eigenfunction corresponding to the eigenvalue Afe of 7, i.e., 70fe — \k4>k- Since the l-pdm is 
selfadjoint, we choose a set of eigenfunctions {0fc}^i which forms an ONB of f). We define the 
operators u : f) -> f) and v : (] — ^ t) by 

U(j)k ■■= {t + lY (l>k ■■= {I + ^kY (t>k and w0fe :=a(l+7)"^(/)fc :=(l + Afc)^5Q,^^^ 

where we abbreviate 1 = 1 f, . We show that 

_ 'w v\ _ ( (1+7)^ a(l+7)"^^ 



V uj \^a*(l+7) 5 (1+7)2 y 

defines a boson Bogoliubov transformation. To this end, we prove the conditions on u and v 
specified in Definition 12.51 We know from the proof of Proposition 15.31 that ([57)) is equivalent to 
([55)) - ([5S)) . Using 07 = 70; and a a* = 7 -|- 7^ we calculate 

uu* — vv* =1+7— (1+7) aa* ~ 1, 

which is the left equation of ([14ap . The right equation of (|14ap is derived from 

u*u-v'^v= (1 +7)^ (1 +7)^ - (1 +7)"'^a^a(l +7)"^ 
and using a^ = a and a a* — -f"^ + j. Furthermore, 

u*v- v^u = (1 +7)5 a{l +7)"^ - (1 -1-7)"^ 0(1+7)2 =0 
because a^ = a and 07 = 70. Thus, we get the left equation of ()14b[) . Analogously, we obtain 

uv^ -vu^ = (1+7)^ (1 +7)^^0^-0(1 +7)"^ (H-7)^ = 0. 
Hence, [/ is a boson Bogoliubov transformation. Since 

tr(w*t;) = tr ('(I + 7)"^ 0*0(1 +7)"^) = tr (7) = tr (7) = w(]N) < 00, 

there is a unitary implementation U : J-^ — > J-^ of the Bogoliubov transformation U by 
Lemma 12.81 

We define a state uju e Z+ by uju{.A) := w(TUAU*) for any A e yl+. We show that the 
generalized l-pdm of wy is 

- - (s :) • 

Since 7(7 = U*^U by Lemma [2.271 ([38)1 is equivalent to 

7 a \ J fO 0\ ^ _ (vv* vu^ 



a* 1 + 7 y V ly \uv* uu 

Thus, we only check that 7 — vv* and o = vu'^ . On the one hand, 

vv* = o (1 + 7) a* = (1 + 7) o o* = 7 
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with a 7 = 7 a and a a* =7 + 7^. On the other hand, 

vu'^ = a(l +7)"^ (1 +7)^ = a. 

Therefore, the Bogohubov transformation U yields Ju — ioi)- I^i particular, we have 

7c/ = (39) 

and juSju = -7(7- 

Next, we show that the only state having jjj as its generalized 1-pdm is the vacuum state. We 
choose {</?n}^i to be a fixed, but arbitrary ONB of f). We consider the set 



K. := Ia' 
and define for any K £ K, 



F U {0} 3no e F, Vn > no : ii:„ = o| 



•'K 



n 

n=l 



{<) 



K„ 



/iO 



and "iiK e J^^ by *^, := a*j^Vl. Note that *(o,o,...) = ^- For any K,L e K with A' ^ L, 
{^K, *l) = 0, {'i'K,^K) = 1, and {'i'KJKeK ^ -^^ forms an ONB of the boson Fock space, the 
so-called occupancy number basis. 

We denote by pu the density matrix corresponding to wc/. With the occupancy number 
representation and the usual Dirac bra-ket notation, we rewrite the density matrix pjj of the 
state LUij as 



Pu 



= Yl Mif.i l*if) (*ih 



K,LeK 



where pk,l '■= i^c/(|^l) {^k\) & C For any n G IN and any K G JC with Kn > 1, we denote the 
vector, in which one of the particles in the state given by (pn is removed, by A" — E^"' , where 
^(") e /C with Et^ = 1 and E^n^ = for all m e IN, m^ n. Then, for all n G F and AT, A e /C 
with Kn > 1, the Cauchy-Schwarz inequality yields 



\^^K,L\ 



^u[ 



|*l)(* 



K-Ei^U 



aa; 



which vanishes since, by p9p . cj[/(a^an) — {'fin, lu Vn) = for all n e F. Consequently, pk,l = 
if any AT, A e /C is different from (0, 0, . . . ), and pKji = 1 for AT = (0, 0, . . . ). As asserted, we 
have 

ojuiA) = {n,An) 

for any A e ^+. Hence, we obtain 

uj{A) ^uju(v*AV) = {vn,Avn) 

for the original state lo. So, w is a pure quasifree state according to Lemma 15.51 D 



30 



We conclude from Lemmas 15.21 and 



Theorem 5.7. Let u: € Z^^n ^^ ^ centered state and 5 = Ij, ® (^ lii) G ^(t) ® ()). Then, the 
following statements are equivalent: 

(i) Lo is a pure quasifree state. 

(ii) The generalized one-particle density matrix 7 0/ the state uj satisfies tr (7) < co and 

jSj = —7. 

Proof. The implication (i) => (ii) is given by the second assertion of Lemma l5.2l and the reverse 
by Lemma 15.61 D 



A consequence of Lemma 15.61 is the following corollary. 
Corollary 5.8. Let lo G Z+ be a state and 



'7 


a 


6' 


a* 


1 + 7 


6 


b* 


r 


1 



the corresponding further generalized 1-pdm with 7 : f) — > f) and a* : t) — )■ [) as defined in Eqs. (J16p 
and (jlSp . respectively. As in (j22l) . the first moment b G \) of the state w is given by {f,b) := 
a;(a(/)) for any f & i). Furthermore, we define the selfadjoint operator Qf : [}©f)®C— >f)©f)©C 
by 



h 





-f 





-111 


f 


-r 


T 


-1, 



for any f €i). Lf 

iQbl^-l, (40) 

then u is a pure quasifree state. 

Proof. If UJ is centered, we have 6 = and (j40p reduces to 7 Qo 7 = ~1 which is equivalent 
to ([57]) . So, Lemma [5.61 directlv yields the assertion. Now, we do not assume the state to be 
centered. Then, for the Weyl transformation Wft : F^ -> J^+, we define the state wq € 2+ 
by wo(A) := uj{WlAWb) for any A £ A+ . First, we show that 60 = 0, 70 = 7 - |&) {b\, and 
ap = a* — \b)(b\ for this state wq. For any / € f), we have 

(60, /) := iOo{a*{f)) = c.(a*(/) - (6, /) 1^) = (6, /) - (6, /) = 0. 

Thus, 60 = and cjq is a centered state. Furthermore, for any /, g G 1], 

(/,7o5) — ^o(a*(.9)a(/)) 

= c.([a*(g)-(fe,g)l^][a(/)-(/,fe)l^]) 
= {f,19)~{f,b){b,g). 
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An analogous calculation yields 

Q-,alg) ■.= uJo{a*{g)a*{f)) 

^oj{[a*{g)~{b,g)t^][a*{f)~{bJ)t^]) 
^u;{a*ig)a*if))^{b,g){bj). 

Next, we consider (l40l) . For every / e f) we decompose the operator Qf as 

Qf^TZfSn}, 

where the operators TZf,S: ()®[}®C— >f)®()®C are given by 

/-Ih 0\ /!(, 

7^/ := ^Ifi and 5:= -If, 

\r T ij Vo -h 

Since TZb is invertible, (1401) is equivalent to 

(7^J77^h)5(7^^77^b) = -7^,*77^6. (gni) 

A straightforward computation yields 

H-\b) {b\ a - \b) (6| _ 0^ 

7^6 7 7^6- a*-|6)(6| li,+7-|6)(5| 

V 1; 

Summarising the results, we obtain 

(70 ao 0^ 

al Ir, +7o 
1, 

which is the further generalized 1-pdm 70 of the state wq. Thus, (PU)) implies 

70^70 = -70- 

Since the upper left 2 x 2- matrix of 70 (which is an operator on [) © t)) is the generalized 1-pdm 
70 and S is diagonal, we find 

70^70 = -70- 

Hence, the generalized 1-pdm 70 fulfills ([57)1 and for the state ljq the requirements of Theorem l5.6l 
are satisfied. Therefore, cjq is a pure quasifree state. D 

An important set of states which are related to the vacuum state via a Weyl transformation 
is the set of coherent states. Recall that a state lo e Z+ is called coherent if there is an / e f) 
and a Weyl transformation W/ : F^ -^ J-^, such that, for any A G A^ , 

uj{A) = {w}n,Aw}n). 

Corollary 5.9. Eq. (I40p is satisfied for every coherent state. 
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Proof. For every coherent state w, we can find G f), such that 

uj{A)^{n,w^Aw;n), 

where W^ : J^"*" — )> T~^ is a Weyl transformation. 
We have b = <j), because, for any f £ i), 

{bj)^{n,w^a*{f)w;n)^{n,[a*{f) + {cpj)t]n), 

where we use (fi, a*{f)ft) = (a(/)il, ft) = 0. For the l-pdm 7 we find 

{f,^g)^{n, [a* (g) + (0, 5) ] [a(/) + (/, ^) l] Q) = (/, 4>) (0, g) = (/, b) {b, g) 
for every /,5 G f) and, thus, 7 = |6) (6|. Furthermore, a* = Ib^lbl by 

(7, a*g) = (r!, [a* (g) + (^, 5) l] [a* (/) + (0, /) l] n) = (0, g) (</>, /) = (/, fe) (6, 5) 
Finally, we obtain 

/O 0> 

nijnb= ih 
Vo 1, 

which obviously fulfills (|^ . D 

5.2. Fermions 

The statements of Sect. lS.ll can also be transferred to fermion systems. The fermion analogue of 
Lemma 15.21 is the following lemma. 



Lemma 5.10. // an operator 7 = ( J. j"-) : () ffi [} — > f) ® f) satisfies < 7 < Ij,, tr(7) < 00, 
and 

f = 7, (41) 

then there is a unique pure quasifree state uj Cz Z^ that has 7 as its generalized one-particle 
density matrix. Furthermore, let uj Cz Z^ be a pure quasifree state. Then, the corresponding 
generalized l-pdm 7 fulfills (|4ip . 

This lemma is a consequence of Theorems 2.3 and 2.6 of [S]. 

Proof. From [5l Theorem 2.3] we conclude that, for every generalized l-pdm 7 , there is a unique 
quasifree state oj E Z^ having 7 as its generalized l-pdm. On the one hand, [5] Theorem 2.6] 
implies that this quasifree state is pure since the corresponding generalized l-pdm is a projection. 
This proves the first assertion of the lemma. 

On the other hand, [SJ Theorem 2.6] also states that the generalized l-pdm of a pure quasifree 
state is a projection which is the second assertion and completes the proof. D 

There is even a one-to-one relation between pure quasifree states and generalized 1-pdms 
fulfilfing (gH). 

Lemma 5.11. Let lu G Z^ . If the generalized l-pdm 7 corresponding to the state uj satisfies 

7' = 7, (42) 

then UJ is a pure quasifree state. 



33 



Proof. Let 1 = If, and 7 = ( J* 1--) be the generalized 1-pdm of w. Eq. (|^ implies a7 = 7a, 
a*7 — 7a*, aa* = 7 — 7^, and a*a = 7 — 7^- We denote by {Aij^j^ the eigenvalues of the 1-pdm 
7 (counting also degeneracies) and choose the corresponding eigenfunctions ^^ G (), z G IN, in such 
a way that {^i}^]^ is an ONB of the one-particle Hilbert space \). Furthermore, let P : t) — > f) 
be the orthogonal projection on the eigenspace of the eigenvalue 1 of 7 and P^ := \jr — P the 
projection orthogonal to P. Note that both projections commute with the 1-pdm and that P is 
also the projection on the eigenspace of the eigenvalue 1 of 7. Furthermore, from 017 = 7a we 
obtain a P () C P [} and a P-^^ C P-'-f). So, P and P-"- commute with 7, 7, a, and a* . 
We define (1-7)^ and (1 - 7)"^ P^ by 

(1 - 7)^ 0. - (1 - A,)^ 4>, and (1 - 7)-^ P^4>, = (1 - A,)"^ P^4>^ 

and consider a Bogohubov transformation U : ()®f)— ?►()©(), C^=(^ir); given by 

m:=(1-7)^ and w := a (1 - 7)"^ P-^ + P. 

First we show that U is indeed a Bogoliubov transformation, i.e., that ([23^, b) hold. The operators 
u and V satisfy 

uu* + vv* = (1 - 7) + [a (1 - 7)"^ pi + pi [(1 _ ^)-5 p^a* + P 



l-7 + a(l-7) ^ P^a* 



P + a (1 - 7) ^ P^P + P (1 - 7) ^ P-i 



7^ , we have 

= 1 - 7 + a a* (1 - 7)"^ P^ + 7P = 1. 



With a*7 = 70;* and a a* = 7 - 
uu* + vv* 

u*u + v'^v — 1 can be shown analogously. Moreover, using (1 — 7)^ P = and (1 — 7)^ P = 0, 

u*w-f u'^u= (1-7)^ a(lL-7)~^P-L + P + {l-'jy^ P^a'^ + P (1-7)^ 

= a P^ + (1 - 7)^ P - P-^a + P (1 - 7)^ 
= 0. 

uv'^ + vu^ = is obtained similarly. Since, furthermore, the (operator valued) entries on the 
diagonal of the matrix U, and, due to a* = —a, those on the off-diagonal as well, are complex 
conjugate to each other, [/ is a fermion Bogoliubov transformation according to Definition 12.361 
The Bogoliubov transformation U has a unitary representation U because 



tr {v*v) 



tr 
tr 



( {1- 7)"^ P^a* +P a (1 - 7)"^ pi 
((1 - 7)"^ 7 (1 - 7) (1 - 7)"^ P^ + P) 



P 



due to a*a = 7 — 7^ and, thus, 

tr {v*v) = tr (7P-L +P) = tr (7) = w(iN) < 00. 

We define a state uju e Z^ by uJui^) '■= uj(V AtJ*) for any A e A^ and denote its density 
matrix by pij. We show that the corresponding generalized 1-pdm 7^/ is given by 



lu = 





1 
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Consequently, 

7r/ = (43) 

and the transformed generalized l-pdm 7[/ is a projection. 

By Lemma [2.551 the Bogoliubov transformation U yields ^u — U*^U. So, U satisfies 



0\ ^ fvv* vu^ 
1 \uv* uu'^ 



that is 7 = vv* and a = vu^ . This is indeed the case since 

vv* ^ P + a{l-jy'^ P^a* =jP + aa*{l- 7)"^ P^ ^ -/P + -/P^ ^ 7, 
vu'^ = a (1 - 7)"^ P^ (1 - 7)5 + a (1 - 7)"^ P^P = aP^ = a. 

Here, we used aP = which follows from Pa*aP = P (7 — 7^) P = 0. 
Let {<Pn}^i denote an arbitrary ONE of f). We define 

/C := |fc e IN ^ {0, 1} 3no e IN Vn > no : fc„ = o| . 

The elements K ^ K, are the occupancy number representations of the fermion Fock space. If 
we define 



-K 



n=l 



for any K ^ K,, the functions ^(0,0,... ) = ^ ^-^d 4'x G -7^^, given by ^^ := c|^51 for A' E K,, K ^ 
(0, 0, . . . ), form an ONB of the fermion Fock space. Furthermore, for any n g IN and any K £ K, 
with Kn = 1, we write K\ {n} for the set where the particle in the state given by (pn is removed, 
but the others are left unchanged. Now, we write the density matrix corresponding to uju as 

K,LeK 

where the coefficients are given by hk,l := ^u{ I^l) (*a'| ) G C for any sets K,L £ K.. Applying 
the Cauchy-Schwarz inequality, we obtain for every n e IN and every pair K,L £ K, with Kn = 1 

\^-Kx\^ = \^u{cl \n) {^\cK\{n}Cn)\ < Wy (cl \^) {il\ CK\{n}CK\{n} |^) (f^| Ci) W[/ (c* C„) . 

By (H5)) . W[/(c*c„) — {(fin, 7c/ ipn) = for every n e IN and fiK.L = if one of the sets K,L £ JC 
is not (0,0,...). Hence, for any A e A^ 

ujuiA) = {n,An) . 

Since [/ is a Bogoliubov transformation with unitary implementation U and invertible, we obtain 
for any A e A" 

uj{A) :^uju{c*Av) = {vn,Avn). 

Therefore, the state lo is pure and quasifree by Remark 12.451 which yields the assertion. D 
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From the last two lemmas we conclude: 
Theorem 5.12. Let lo Cz Z^ be a state. Then the following statements are equivalent: 
(i) u) is a pure quasifree state, 
(ii) The generalized 1-pdm 7 of the state uj satisfies 

7^ =7. 

Proof. The implication (i) => (ii) is given by the second assertion of Lemma lS.lOl and the reverse 
by Lemma [5. Ill D 

A. Generalized 2-pdm as 7 x 7-Matrix 

In this appendix, we give a more explicit, basis dependent form of the generalized 2-particle 
density matrix F. We assume {(f>k}'^^i to be a fixed, but arbitrary ONB of f). Recall that F is 
defined as a 7 x 7-matrix on [}sim in Definition l2.32l In order to simplify notation, it is convenient 
to define some operators and functionals. 

Definition A.l. Let /i, /2,5i,52, /,5 e f], and ^ e C. We define 7^(3) := {F e f) (g) 
i)\ J2tLi {^k ® 4>k^F)^^^^ < 00} C [) (g) () and the following linear maps: 

H , ((7i(g52,Ai(/i(g)/2)) ■.= u{a*{fi)a*{f2)a*{g^)a{gi)), 
§fl , (51 ® 52, a; (/i 0/2)) :-^(a*(/i)a*(/2)a*(52K(gi)), 
[} , {gi® 5-2, A (/i ® /2)) := i^ {a*(fi)a*(-g^)a{g2)a(f^)) , 

(g, A, (/i (g> /2)) := LU (a*(/i)a*(/2)a(g)) , 

(.9,^2 (/i ® /2)> := io {a*{f,)a*{h)a*m , 

{9,Qi (/i ® /2)> := c. (a*(/i)a(72)a(5)) , 

{9,Q2 (/i ® /2)> := LO (a*(/i)a(72)a*(5)) , 

00 

>l} , B -.^ Y^ \(|3^ (® 4'^) {(/Jk (® (l^kl , 
i,k=l 

00 

/3i: I?(B)^C, /3i:=^(0,®0,|. 



Ai 


[}(Xit) ^ f) 


A^ 


f)(g)fH- f) 


A: 


()«)() ^ [}( 


Ai 


t)(X) f) ^ f) 


A*2 


f](Xif) ^ f} 


Qi 


f)® f) ^ t) 


Q2 


[}® f) -^ f) 


B : 


p(i?) ^ () 


P2: 


C(8)() ^ f) 



Furthermore, recall from Remark 12.291 that h is given by (6, /)t := w {a*{f)). 

As we already pointed out in Proposition 12.331 the generalized 2-pdm is selfadjoint. Since 
therefore F^ = Tji, it suffices to state the entries F.^ for 1 < i < j < 7. Using the notation 
specified before, we have: 

fii =F : ()(g)f) -> f)(8)[), 

fi2 = A* : f) (g) () ^ f] (g) [}, 

fi3 = A*Ex +{a®t)B : f)(g[)^f)(gf). 
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ri4 = 


= A2: i)«)[)^fl(g)(), 


ri5 = 


= Al: ()-^i)®l), 


fl6 = 


^ A2 : [) -^- f) (g) f), 


fl7 = 


= {a <g) 1) l3l : C -^ f)(g)f), 


T22-- 


= ExA + 7«)1: i)<g)\)^i)<g)i), 


r23 = 


= A* + (7 ® 1) (S + Ex) : () ® [) ^ [) ® (), 


f24 = 


= ExAi + (a «) 1) (1 + Ex) : () «) f) ^- i) (g) (), 


f25 = 


- Qt : f) ^ d ® f), 


r26 = 


-Q;: f)^f)(gf), 


r27 = 


= (l(g7)/3* : C^(i(gf), 


Tsa = 


== AEx +(l(g)7)B + B(l(g)7) + B + l®7: [)®[)^f)(8)[) 


f34 = 


= Ai + (1 ® a) (1 + Ex) + B (1 (g) a) : I) I) ^ f) f), 


f35 = 


^ExQl+Plb* : f)^f)®f), 


Tae = 


= ExQ;+/3*r : f)^f)®f), 


r37 = 


= (1 (g) 1 + 1 (X) 7) /?! : C ^ () «) f), 


r44 : 


= r^ + (l(8)l+7(8)H-lL(g)7)(l + Ex):[)(8)f)-^f)(g)f), 


r45 = 


= 32 : f) -> f) «) f), 


r46 = 


= 3i + (1 + Ex) (1 ® fe) /32* : f) -> f) f), 


r47 = 


= (a*(Xl)/3* : C^f)(8)f), 


r55 = 


= 7: f) ^ t), 


r56 = 


= a : ft ^ f], 


r57 = 


= 5: C^(), 


Tee = 


= 1 + 7: [)^tl, 


r67 = 


= 5: C^t), 


r77 = 


= 1: C^C. 



Remark A. 2. Both generalizations of the l-pdm, given in Definitions l2.24l and l2.281 are contained 
in the generahzed 2-pdm, namely 





r \ 






r,5 


Tfis 


r75 


55 


and 


7 = 


r.fi 


Tee 


r76 


56 






[h7 


r67 


r77 
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